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c^J ' Context. We present the discovery of very high energy (VHE, E > 100 GeV) y-ray emission from the BL Lac object 1ES 1215+303 by the MAGIC 
^* , telescopes and simultaneous multi-wavelength data in a broad energy range from radio to y-rays. 
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Aims. We study the VHE y-ray emission from 1ES 1215+303 and its relation to the emissions in other wavelengths. 

Methods. Triggered by an optical outburst, MAGIC observed the source in 2011 January-February for 20.3 hrs. The target was monitored in the 
optical R-band by the KVA telescope that also performed optical polarization measurements. We triggered target of opportunity observations with 
' . . , the Swift satellite and obtained simultaneous and quasi-simultaneous data from the Fermi Large Area Telescope and from the Metsahovi radio 
^» . telescope. We also present the analysis of older MAGIC data taken in 2010. 

' Results. The MAGIC observations of 1ES 1215+303 carried out in 2011 January-February resulted in the first detection of the source at VHE 
1 with a statistical significance of 9.4 cr. Simultaneously, the source was observed in a high optical and X-ray state. In 2010 the source was observed 
, in a lower state in optical, X-ray, and VHE, while the GeV y-ray flux and the radio flux were comparable in 2010 and 201 1. The spectral energy 
distribution obtained with the 201 1 data can be modeled with a simple one zone SSC model, but it requires extreme values for the Doppler factor 
or the electron energy distribution. 

Key words. Gamma rays:galaxies-BL Lacertae objects:individual: 1ES 1215+303 



1. Introduction 

Most of the extragalactic sources from which very high en- 
ergy (VHE, > 100 GeV) y-ray emissions has been detected are 
blazars. These objects are commonly believed to be a subtype 
of active galactic nuclei (AGN) whose relativistic jet points very 
close to the line of sight of the observer. Blazars are character- 
ized by high amplitude variability at all wavebands from radio 



to y-rays. The correlations between the different energy bands 
are complicated, but in general it seems that high states in lower 
energy bands (e.g. optical) are accompanied by high states in the 
higher energies (i.e. y-rays) at least in some sources (see e.g. 
lReinthaletalll20L2l) . 

The spectral energy distribution (SED) of blazars exhibits a 
generic two-bump structure: one peak with a maximum in the 
spectral range from radio to X-rays and a second peak in the in- 
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terval from X-ray to y-ray. The radiation is produced in a highly 
beamed plasma jet and the double peaked SED is often explained 
by a single population of relativistic electrons. The first peak is 
due to synchrotron emission in the magnetic field of the jet and 
the second peak is caus ed by invers e Compton (IC) scattering 
of low-energy photons dReesl 1 19671) . The low-energy photons 
can originate externally to the je t (external Compton scatter- 
ing, |Dermer^Schlickeise3, [l993|) or be produced within the jet 
via s ynchrotron radiation (s ynchrotron self-Compton scattering, 
SSC. lMaraschi et allll992l) . 

Blazar is a common term used for Flat Spectrum Radio 
Quasars (FSRQs) and BL Lac objects (BL Lacs), which are 
thought to be intrinsically different. The FSRQs show broad 
emission lines in their optical spectra while the BL Lacs have 
featureless spectra with weak or no emission lines possibly 
masked by a strong emission from the jet. This indicates 
that in BL Lac objects the main population of seed photons 
for Compton scattering should originate from the synchrotron 
emission. Indeed most of the SEDs of BL Lacs are well de- 
scribed with simple SS C model (e.g. lBloom &Marscherlll996t 
iTavecchio et~ailll998l) . 

MAGIC has been successfully performing optically trig- 
gered VHE y-ray observations of AGN since the start of its 
operations. The triggers have been provided by the Tuorla 
blazar monitoring program^ and the target of opportunity (ToO) 
observations with MAGIC have resulted in the discovery of 
five new VHE y-ray emit ting sources (Mrk 180, lAlbert et all 
(I2006bl): 1ES 1011+4 96, lAlbert et al.l d2007l): S5 0716+714 , 
lAnderhub et ail d2009h ; B3 2247+381, lAleksic et al.l d2012al) : 
and 1ES 1215+303, this paper). However, in many cases it has 
not been possible to confirm with high statistical significance if 
the sources were in higher VHE y-ray state than usual during 
the observations. The lo ng-term studies of individual VHE y- 
ray blazars like Mrk 4 2 1 (lAcciari et all 1201 ll) and PG 1553+113 
( lAleksic et al.Ll2012bl) have also yielded controversial results on 
the correlation between the two broad energy ranges. Thus, to 
date, the connection between the optical and VHE y-ray states 
has remained an open question. 

1ES 1215+303 (also known as ON 325) is a high syn- 
chrotron peaki ng BL Lac obje ct ( Abdo et all 1201 Obi) with red- 
shift z = 0.130 (|S^amaetail|2^|niowewr,z = 0.237 isalso 
reported in the literature, e.g. NEE0). The source was classified 
as a promising candidate TeV blazar dCostamante & Ghisellinil 
120021 ITavecchio etail 120101) and has been observed several 
times in VHE y-rays before the observations presented here. 
The previous observations yielded only up per limits, Whipple : 
F(> 430GeV) < 1.89 x 10 11 cm 2 s"\ (Koran et all 120041) : 
MAGI C: F(> 120GeV) < 3.5 x 10 11 cm 1 s~\ ( lAleksic et all 
2011a). The source was list ed in the Fermi Lar ge Area Telescope 
(LAT) bright AGN catalog dAbdo et alll2009l) as showing a hard 
spectrum (F = 1.89 ± 0.06). It underwent a large outburst in 
late 2008, and in this catalog 1ES 1215+303 is the only high 
energy peaking source that show s significant variability. In the 
second Fermi-L AT AGN catalog dAckermann et all 1201 lb other 
high synchrotron peaking sources have also been flagged as vari- 
able. 

In the first days of 2011 January 1ES 1215+303 was ob- 
served to be in a high optical state. This triggered MAGIC ob- 
servations, extending until 2011 February, that resu l ted in the 
discovery of VHE y-rays from the source dMariottiL l20"llt) . In 
this paper we present the results of the 2011 January-February 



1 http://users.utu.fi/kani/lm/ 

2 http://nedwww.ipac.caltech.edu/ 



observations. We also present the previous observations of 
1ES 1215+303 with the MAGIC telescopes performed in 2010 
January-February and 2010 May- June that produced only a hint 
of signal. For all epochs we present simultaneous and quasi- 
simultaneous multi-wavelength data from radio, optical, X-ray, 
and GeV y-rays. 



2. Observations and Data Analysis 

The observations of 1ES 1215+303 were performed in a broad 
wavelength range (from radio to VHE y-rays) by 5 different in- 
struments. This is the first time that such a broad wavelength 
range is covered for this source in quasi-simultaneous observa- 
tions. 



2.1. MAGIC 

MAGIC consists of two 17 m Imaging Air Cherenkov 
Telescopes (IACTs) sensitive to y-rays with energy above 
50 GeV in standard trigger mode (which is the lowest trigger en- 
ergy threshold among the existing IACTs). The system is located 
in the Canary Island of La Palma, 2200 m above sea level. Since 
fall 2009 the telescopes are working together in stereoscopic 
mode which ensures a sensitivity of < 0.8% of th e Crab Nebula 
flux ab ove 300 GeV in 50 hrs of observations dAleksic et all 
The field of view of the each MAGIC camera has a di- 
ameter of 3.5°. 

1ES 1215+303 was observed by MAGIC in 2010 January- 
February, 2010 May -June , and 201 1 January-February, for a to- 
tal of 48 hrs. The observations were done in the so-called wobble 
mode (i.e. with the source offset by 0.4° from the camera cen- 
ter), which provides a si multaneous estimat e of the background 
from the same data set dFominet all 11994 . While most of the 
data were taken in dark night conditions, a small fraction were 
taken in presence of moderate moonlight. The data span a range 
of zenith angle from 1 ° to 40° with most of the data taken be- 
low 25° (in 2010 the mean zenith angle was ~ 13°, and in 201 1 
~ 8°). 

The data were anal yzed using the standar d MAGIC software 
and analysis package dAleksic et all l2012d). A nother VHE y- 
ray emitter, 1ES 1218+304 (lAlbert et alll2006al) is present in the 
same field of view as 1ES 1215+303. The sources are separated 
by ~ 0.8°, which is much larger than the point spread function 
(PSF) of the MAGIC telescopes (~ 0.1°), so there was no source 
confusion or contamination. However, these sources have nearly 
the same Right Ascension, so in the standard wobble setup used 
in the 2010 January-February observations, the background esti- 
mation region partially overlapped with the 1ES 12 18+304 posi- 
tion. This would result in an overestimate of the background, so 
this region was excluded from the background estimate. In the 
later observations (2010 May- June, and 201 1 January-February) 
the wobbling offset direction was changed to have the standard 
background estimation regions far from the second source. 

After the data quality selection, based mainly on the rate of 
stereo events, the data samples of January-February 2010, May- 
June 2010 and January-February 2011 contain 19.4, 3.5, and 
20.6 hrs of good quality data respectively. Because of the differ- 
ent positions of the source in the camera, and the variable nature 
of AGN, we decided to split the analysis into these 3 periods. 
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2.2. Fermi-LAT 

The Fermz-LAT is a pair conversion telescope designed to cover 
the energy band from 20 Me V to greater than 300 GeV. It oper- 
ates in all-sky survey mode and therefore can provide observa- 
tions of 1ES 1215+303 simultaneous to MAGIC. In this paper 
the standard LAT Science Tools (version v9r23p0) were used to 
analyze the data collected in the time interval from 2008 August 
5 to 201 1 March 22. For this analysis, only events belonging to 
the "Diffuse" class (which have the highest probability of being 
photons) and located in a circular Region Of Interest (ROI) of 7° 
radius, centered at the position of 1ES 1215+303, were selected 
(using Pass 6 event selection). In addition, we applied a cut on 
the zenith angle (< 100°) limb y-rays and a cut on the rocking 
angle (> 52°) to limit Earth limb contamination. 

The data analysis of 1ES 1215+303 is very challenging 
due to the presence of several y-ray emitting sources in the 
same ROI. 1ES 1218+304 is located at a distance of just 0.8° 
from the source of interest. Another well known VHE emit- 
ter, W Comae, is located at ~ 2° from the latter source. 
Thus, the LAT analysis was restricted to energies above 1 GeV 
where the Fermi-LAT PSF is sufficiently narrowQ to separate 
1ES 1215+303 from the other sources in the ROI. The un- 
binned likelihood method was applied to events in the en- 
ergy range from 1 GeV to 300 GeV. All point sources from 
the 2FGL catalog (iNolanetall l2012h located within 12° of 
1ES 1215+303 were included in the model of the region. 
Sources located within a 5° radius centered on 1ES 1215+303 
position had their flux and photon index left as free parame- 
ters. The diffuse Galactic and isotropic components (including 
residual instrumental background) were modeled with the pub- 
licly available files gll_iem_v<82_P6_Vll_DIFFUSE. fit and 
isotropic iem_v(82_P6_Vll_DIFFUSE.txtQ. The normalizations 
of the components comprising the total background model 
were allowed to vary freely during the spectral point fitting. 
The instrument response functions P6_V1 1_DIFFUSE were used. 
The successful separation of flux between 1ES 1215+303 and 
1ES 1218+304 was verified by the absence of any significant 
correlation between their light curves. The systematic uncer- 
tainty in the flux is estimated as 5% at 560 Me V and 20% at 
lOGeVandaboveB 

2.3. Swift 

The Swift satellite dGehrels e t al., 2004j) is eq uipped with three 
telesc opes, the Burst Alert Telescope (BAT; iBarthelmv et al.L 
l200l coveri ng the 15-150keV e nergy range, the X-ray tele- 
scope (XRT; iBurrows et al.L 120051) covering the 0.2-10keV en- 
ergy b and, and the UV/Optical Telescope (UVOT; lRoming et al.L 
120051) covering the 1 80-600 nm wavelength range. 

A Swift ToO request was submitted on 2011 January 3. 
The Swift observations started on January 4 until January 12 
with four ~ 5 ks exposures in photon counting mode. The data 
were processed with standard procedures using the FTOOLS 
task XRTPIPELINE (version 0.12.6) distributed by HEASARC 
within the HEASoft package (v.6.1 0). Events with grade s 0- 
12 were selected for the data (see IBurrows et all l2005l) and 
the latest response matrices available in the Swift CALDB 
(v.20 100802) were used. For the spectral analysis, we extracted 



the source events in the 0.3-10keV range within a circle with a 
radius of 20 pixels (~ 47 arcsec). The background was extracted 
from an off-source circular region with a radius of 40 pixels. 

The spectra were extracted from the corresponding event 
files and binned using GRPPHA to ensure a minimum of 25 
counts per energy bin, in order to obtain reliable x 1 statistics. 
Spectral analysis was performed using XSPEC version 12.6.0. 
The neutral hydrogen-equivalent column density was fixed to 
the Galactic value in t he direction of the source 1.74 • 10 20 cm -2 
dKalberla et al.L 12001 . 

Swift/UVGT observed the source with all filters (V, B, U, 
UVW1, UVM2, UVW2) for four nights. UVOT source counts 
were extracted from a circular region 5 arcsec-sized centered on 
the source position, while the background was extracted from a 
larger circular nearby source-free region. These data were pro- 
cessed with the uvotmaghist task of the HEASOFT package. 
The observed magnitudes h ave been corrected fo r Galactic ex- 
tinction E B - y = 0.024 mag dSchlegel et all 1 1998b . applying the 
for mulae by Pejll 1992h and finally converted into fluxes follow - 
al|<12 



ing lPoole et all (120081) 



2.4. KVA 

The KVA optical telescopes are located in La Palma, but are 
operated remotely from Finland. The two telescopes are attached 
to the same fork. The larger telescope has a mirror diameter of 
60 cm and the smaller 35 cm. 

The 35 cm telescope is used for simultaneous photomet- 
ric observations with MAGIC, but also to monitor potential 
VHE y-ray candidate AGN in order to trigger MAGIC obser- 
vations if the sources are in high optical states. The observa- 
tions are performed in the R-band and the magnitude of the 
source is measured from CCD images using differential pho- 
tometry. The comparison star m agnitudes for 1ES 1215+303 
are from iFiorucci & Tostil (Il996l) . and the mag nitudes are con - 
verted to flux using the formula and values from lBesselll d 1979b . 
1ES 1215+303 has been observed regularly as part of the Tuorla 
blazar monitoring program since 2002. 

The 60 cm telescope i s used for polarime t ric obs erva- 
tions (see e.g. iPiirola et all I2005L lAleksic et all 1201 lbl) . For 
1ES 1215+303 polarimetric observations were performed on six 
nights from 2011 January 7 to January 17. The degree of po- 
larization and position angle were calculated from the intensity 
ratios of the ordinary and extraordinary beams using standard 
formulae and semiautomatic software specially developed for 
polarization monitoring purposes. 

2.5. Metsahovi radio telescope 

37 GHz radio observations were made with the 13.7m 
Metsahovi radio telescope located in Kylmala, Finland. The tele- 
scope, the observation methods, and the data analysis proce- 
dure are described in e.g. lTerasranta et aTl (Il998l) . The telescope 
detection limit is ~ 0.2 Jy under optimal conditions and since 
1ES 1215+303 is a rather weak source at 37 GHz it can only be 
observed under good weather conditions. Typically, an accept- 
able measurement of the source is obtained approximately once 

3 http://www.slac.stanford.edu/exp/glast/groups/canda/archive/pass6vl iper month. Data were obtained simultaneously with the MAGIC 
/lat_Performance.htm observations in 2010 June, but in 2011 January-February the 

4 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.htmlweather did not allow simultaneous observations with MAGIC, 

5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html the closest points being from 2010 December and 201 1 March. 



3 



J. Aleksic et al.: 1ES 1215+303 



3. Results 

3.1. MAGIC results 

The MAGIC data were divided in three samples correspond- 
ing to three observation epochs: 2010 January-February, 2010 
May-June, and 2011 January-February. The so-called 8 2 plots 
(the distribution of the squared angular distance between the ar- 
rival direction of the events and the real position), for energies 
above 300 Ge V, corresponding to the three observation epochs, 
are shown in Fig. Q] The computation of the number of the ON 
(signal) and OFF (background) events was performed in a fidu- 
cial signal region of 6 2 < 0.01 deg 2 , and using 5 background 
regions (4 in case of the January-February 2010 data). In 2010 
January-February (left panel) 194 ON events were detected over 
144.6 + 6.0 OFF events, wit h a significance level of 3.5 cr (us- 
ing Eq. 17 in lLi & Malll983l) . In 2010 May-June (middle panel) 
the observation time was much shorter and no excess events 
were present. For 201 1 data (right panel) the numbers were 251 
ON over 119+ 4.8 OFF corresponding to a ~9.4 cr significance, 
which is the first detection of VHE y-rays from this source. 

In Fig. |2 we show the significance map of the sky region 
for energies above 300 GeV for the 2010 (January-February and 
May-June combined) and 2011 observations. 1ES 1218+304 is 
clearly visible in both maps while 1ES 1215+303 was fainter in 
2010 than in 201 1. The 1ES 1218+304 data analysis and results 
will be addressed in a separate paper. 

For the light curve and spectrum determinations softer cuts 
were applied that have a higher y-ray efficiency. The light curve 
(in a 5 -days bins) above 200 GeV of the 2011 data is well de- 
scribed by a constant flux of (7.7 + 0.9) x 10~ 12 cnT 2 s _1 Ct 2 /«dof 
= 0.56 / 3), which corresponds to about 3.5% of the Crab Nebula 
flux. Assuming that the hint of a signal seen in the 2010 data 
is a y-ray excess the corresponding flux was F(> 200GeV)= 
(3.4 ± 1.0) x 10 _12 cm _2 s _1 , which is less than half of the flux 
measured in 20 1 1 . The hypothesis of constant flux between 20 1 
and 201 1 is excluded at the level of 3.1 cr. 

The derived VHE y-ray spectrum for the 201 1 observations 
can be described by a single power law (^ 2 /nd f = 5.2/3, see 
Fig.0: 
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in the fitting range 70 GeV - 1.8 TeV Since the spectral in- 
dex of 1ES 1215+303 is similar to that of the Crab Nebula 
and the source is relatively bright, we can direct ly use the 
systematic errors estimated in Aleksic et al] d2012d) . The sys- 
tematic error of the slope is +0.15 and in the energy range 
of the 1ES 1215+303 spectrum, the error in the flux normal- 
ization (without the energy scale uncertainty) was estimated 
to be 11%. The systematic error in the energy scale is 15%. 



EBL models dDoi 


-mnsuez et al. . 201 1 : Kneiske & DoleL 1201(1 


Franceschini et al. 


, 2008; Primacketal. 


12005b and the maxi- 


mum high UV EBL model described in 


Albert et al.l (120081) for 



z = 0.130. The results are shown in Fig. [3] As denoted in the 
Figure by the shaded area , at this redshift th e EBL models agree 
well. The EBL model of iDommguez et aTl d201 ll) was used to 
calculate the final intrinsic spectrum since this model is based 
on an observational approach. 

3.2. Fermi-LAT results 

The light curve of 1ES 1215+303 was obtained in the energy 
range from 1 to 100 GeV, in 14-day bins from 2008 August 



to 2011 March (F ig. 0J. It shows the major flare reported in 
lAbdo et all d2009l) at the beginning of the Fermi mission. There 
is a hint of enhanced flux during 2010 November (MJD 55500, 
duration only one bin, i.e. 14 days) but very little variability oth- 
erwise, especially at the two MAGIC observation epochs (2010 
January-June and 2011 January-February). To maximize the 
number of photons the spectral energy distribution was derived 
using the whole MAGIC dataset (2010 January- June and 2011 
January-February). The spectral energy distributions are shown 
in Fig. In 2010 January- June the integral flux, F (1 - 100 GeV) 
is (4.9 + 0.7) x 10~ 9 cirrV and the photon index 2.1 + 0.1, 
while in 2011 January-February F(l - 100GeV)= (7.3 + 1.6) x 
10~ 9 cm" 2 s _l and the photon index 2.0 + 0.2. The mean detected 
flux was ~ 50% higher in 2011 January-February than in 2010 
January-June, but due to large error bars the increase was not 
statistically significant. The spectral index was constant within 
the error bars. 



3.3. Swift results 

The results of the Swift/XKT observations are summarized in 
Table 1. The source showed the highest flux on 2011 January 
8 (MJD 55569.1) and previous/subsequent observations from 
2009 December (MJD 55168.7)/ 2011 April (MJD 55674.2) 
show significantly lower flux. For the X-ray spectra both log 
parabola (in the form ~ £- a - M °g( E \ with E being the energy 
in keV) and a simple power-law fit were tested. The best fit was 
achieved with a log parabola law model in the range 0.3-10 keV 
for four observations while a simple power law, in the range 
0.5-10 keV, provided a better fit for three of the observations. 
Generally, a log parabolic fit suggests that there is curvature in 
the X-ray spectra but for 1ES 1215+303, the difference between 
log parabolic and power law fits is small so no strong conclu- 
sions can be drawn. Because of the different fits a comparison 
between the spectral slopes is difficult, but for the highest flux 
night the spectral index is marginally harder than for the low 
state observations. 

The Swift/WOT results from 2011 January ToO observa- 
tions show constant brightness with V-band magnitude = 15.06+ 
0.10, B= 15.38 + 0.10, U= 14.53 + 0.08, UVW1 = 14.43 + 0.08, 
UVM2= 14.35 + 0.06, and UVW2= 14.46 + 0.06. However, in 
all bands the source is clearly brighter than in the previous ob- 
servation (2009 December: V= 15.60 + 0.10, B= 15.95 + 0.10, 
U= 15.12 + 0.08, UVW1= 15.07 + 0.08, UVM2= 15.00 + 0.06, 
and UVW2= 15.15 ±0.06). 



3.4. KVA and Metsahovi results 

In the optical R-band the source is clearly variable on daily 
and yearly time-s cales. The host gala xy contributes a flux of 
0.99 + 0.09 mJy dNilsson et al.L 120071) and when this contribu- 
tion was subtracted from the measured flux, the AGN core was 
found to be ~ 40% brighter in 2011 January-February (average 
total flux 3.64 mJy) than in 2010 January-June (average total flux 
2.55 mJy). Similarly, it was found that during the 201 1 January- 
February observations the flux varied by ~25% (core flux be- 
tween 3.2 mJy and 4.1 mJy). 

During 2011 January the optical polarization was ~ 9% 
while during the follow up observation in 2011 April it was 
higher, ~ 15%. The position angle (PA) was only slightly vari- 
able between 140 and 150 degrees. 

In the radio band the source is rather weak and does not 
show strong variability. The 37 GHz flux from the Metsahovi 
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Fig. 1. Distributions of the (r parameter for 1ES 1215+303 signal (black histograms) and background estimation (gray histograms) 
for the three observation periods: January-February 2010 (left), May-June 2010 (middle), and January-February 201 1 (right). The 
vertical dashed line corresponds to the apriori defined signal region 2 < 0.01 deg 2 . 



MJD 


Obs. Time [ks] 


Flux (2-10keV) [10^'- erg/cm 2 /s] 


a (T for PL) 


b 




55168.6799 


4.99 


1.21+0.19 


2.56 + 0.10 


0.34 ± 0.34 


1.19/25 


55565.0340 


4.39 


2.74 ± 0.25 


2.41 ± 0.08 


0.37 ± 0.24 


0.87/42 


55569.1281 


2.38 


3.02 ± 0.40 


2.29 ± 0.16 




1.23/18 (PL) 


55571.1327 


4.07 


1.69 + 0.17 


2.65 ± 0.14 




1.23/18 (PL) 


55572.1361 


4.27 


1.45 ± 0.20 


2.64 ± 0.09 


0.28 ± 0.27 


1.15/32 


55573.1396 


2.99 


1.73 ± 0.25 


2.46 + 0.11 


0.66 ± 0.37 


1.26/25 


55674.2438 


2.34 


1.30 + 0.30 


2.67 ± 0.25 




0.48/8 (PL) 



Table 1. Data summary and results for the Swift/XRT ToO observations. The datasets in the first/last rows are prior/subsequent to 
the MAGIC observations and are reported for comparison. For each dataset the following quantities are reported: the MJD time 
of the beginning of the observations; the exposure time; the integral flux in the 2-10 keV band; the a and b parameters for the log 
parabola fit (or the photon index F in case a simple power-law is used, see text); the reduced^- 2 with number of degrees of freedom 
Wdof • PL indicates when the simple power law is used instead of the log parabola. 



radio telescope has a similar level (0.3-0.4 Jy) in 2010 and 201 1, 
although there were no radio observations during 201 1 January- 
February. 



4. Interpretation 

In this section we discuss the quasi-simultaneous light curves, 
showing how they establish connections between different en- 
ergy regimes and locate the emission region. The spectral energy 
distribution is reconstructed for the first time from radio frequen- 
cies to TeV energies for 1ES 1215+303, allowing us to study the 
capability of the one-zone synchrotron self Compton model to 
reproduce the constructed SED. 

4.1. Multi-wavelength behavior 

The long-term multi-wavelength light curve, from radio to VHE 
y-rays, is shown in Fig. [6] The MAGIC light curve shows a lower 
flux in 2010 (January-February and May- June) than in 201 1 (by 
a factor of 2). The large uncertainties in the Fermi-LAI mea- 
surement do not allow us to conclusively say whether a simi- 
lar flux enhancement also occurred in the 1 - lOOGeV energy 
range (see section 3.2). In X-rays the source was in a high state 
(enhanced by a factor of 2) in 201 1 January compared with pre- 
vious observations. In the optical band the average flux during 
MAGIC observations in 2010 was 3.5 mJy, while in 201 1 it was 
4.6 mJy. Thus, the source was clearly in outburst during early 
201 1, at least in VHE y-rays, X-rays, and the optical band. There 
were no simultaneous radio observations, but both the previous 
and subsequent observations showed low flux, suggesting that 
the outburst might have originated rather close to the central en- 
gine where the emission region is opaque at radio wavelengths. 



However, as the simultaneous observations are missing the exis- 
tence of a simultaneous radio flare cannot be excluded. 

During the 201 1 January-February observations (Pig. [7j, the 
MAGIC light curve is consistent with a constant flux. The source 
was in a rather low state in the Fermi-LAT energy range and no 
short term variability was detected. In X-rays and optical the 
source was variable during the MAGIC observations: the first 
two X-ray exposures gave a higher flux than for the latter three. 
The X-ray spectra show hints of hardening with higher flux, but 
they are statistically the same. The MAGIC observations started 
when the optical flux was decreasing, but during 2011 January 
the optical light curve showed several small flares. The X-ray 
light curve was more sparse and showed only one flare, but the 
comparison of simultaneous optical and X-ray points shows the 
same pattern in the light curves, indicating that the X-ray and 
optical emissions originates from the same region. 

In addition to multi-wavelength variability studies, the op- 
tical polarization measurements have proven to be a powerful 
tool to analyze the e mission scenarios in the blazar jets (e.g. 
iMarscheret all 120081) . Polarization traces the magnetic field of 
the jet. A net polarization oriented either parallel or perpendic- 
ular to the projected jet axis can be confused by shocks and the 
signatures are visible in optical polarization. The optical polar- 
ization measurements from 2011 January show little variabil- 
ity in polarization degree (average ~ 9%) or PA (varying be- 
tween ~ 140° - 150°) during the MAGIC observations, but the 
follow-up observations from 2011 April (Pig. [6]) show a higher 
polarization, ~ 15%. Unfortunately, the polarization observa- 
tions missed the peak of the first optical outburst and our data 
sample is very small. Ilkeiiri et all (1201 ll) monitored the photo- 
polarimetric behavior of the source in 2008-2009 and their ob- 
servations seem to show similar polarization trends (i.e. a de- 
creasing polarization during outbursts). They also found that the 
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Fig. 2. Significance maps (> 300 GeV) from MAGIC observa- 
tions performed during 2010 January-February and 2010 May- 
June (combined together, total time 22.0 hrs, top) and 2011 
January-February (total time 20.3 hrs, bottom). 
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Fig. 4. Long-term light curve of 1 ES 1 2 1 5 +303 from Fermz-LAT 
between 1 and 100 GeV. 
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Fig. 5. Spectral energy distribution from Fermi-LAH derived for 
2011 January-June (top) and 2011 January-February (bottom). 
The upper limits have be en computed when the test statistics (see 
e.g. iMattox et al.L 1 19961) in the energy band were lower than 4. 
The bow-ties are derived from the unbinned likelihood analysis. 



Fig. 3. Observed and deabsorped VHE y-ray spectra for a red- 
shift of 0.130. The EBL model of iDommguez et al l (1201 lh was 
used, the gray area shows the spread of the EBL models. The 
arrow shows the systematic error of the measurement. 



PA was almost constant at ~ 150°, which agrees wit h our obser- 
vation s and with the historical data from 1981-1989 dWills et al.L 
1201 lh showing PA values from ~ 130 - 170°. Such preferred 
posit ion angles have bee n observed for several BL Lac objects 
(e.g. lJannuzi et alj. 1 19941) and implies long-term stability of the 
structure of the region producing the polarized emission e.g. the 



existence of a optical polarization core. In first order, if the opti- 
cal outburst was produced by a shock traveling along the jet, one 
would expect the polarization degree to increase during the out- 
burst. However, if there is a standing shock (optical polarization 
core) present, another shock with a different magnetic field ori- 
entation colliding with the standing component could produce an 
outburst in the total flux, but de crease the observed level of po- 
larization (IVillforth et al.L l2010l) . A detailed photo-polarimetric 
study based on more data would be needed to further test this 
hypothesis. 
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Fig. 6. Long-term multi-wavelength light curve of 
1ES 1215+303 from 2009 December to 2011 May. The 
vertical line shows the beginning of the MAGIC 2011 obser- 
vation campaign, a) In the MAGIC light curve 2011 data are 
binned in 5-day intervals. 2010 data are divided in January- 
February and May-June bins, b) The Fermi-LAT light curve 
(1 - lOOGeV) has bins of 14 days and the points with arrows are 
upper limits, c) The SwiftPiKT light curve is derived from the 
target of opportunity observations performed during the MAGIC 
observations and archival data, d) The R-band light curve shows 
hourly average flux of the source, the error bars are smaller than 
the symbols in most cases, e) The optical polarization (filled 
circles, left axis) and polarization position angle (triangles, right 
axis) are hourly averages, f) 37 GHz radio light curve from the 
Metsahovi radio observatory. 



4.2. Spectral energy distribution 

The SED of 1ES 1215+303 in both MAGIC observation epochs 
is shown in Fig. 8. The 201 1 high energy bump is constructed us- 
ing the MAGIC deabsor bed spectrum (using the EBL model of 
iDommguez et al.L 1201 lb and the simultaneous Fermi-LAT spec- 
trum (collecting all photons from 2011 January-February). As 
stated in section 4.1, the low energy bump was variable dur- 
ing the period and is constructed for the night MJD 55569 
that showed the highest Swift flux and for which there are si- 
multaneous KVA and UVOT observations. The contribution of 
the host galaxy was subtracted from the R-band flux following 
iNilsson et al.l (l2007h . The host galaxy also contaminates the V, 
B, and U bands of the UVOT data, but its contribution should be 
negligible in the UV. As we have no direct measurements of the 
host galaxy contribution in V, B and U bands we extrapolated 
the magnitudes from the R-ba nd value using t he ga laxy colors 
of elliptical galaxies at z = 0.2 dFukugita et al l fT995h . 

For the 2010 MAGIC data set, we could not derive a spec- 
trum because of the low significance of the signal but we report 
the flux between 300 GeV and 1 TeV (assuming the same spec- 
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Fig. 7. Multi-wavelength light curve of 1ES 1215+303 from 
2011 January to February, a) In the MAGIC light curve, the 
data are binned in 5-day intervals, b) The Fermi-LAT light curve 
(1 - 100 GeV) has bins of 14 days, c) Swift/XRT light curve, d) 
UVOT optical and UV light curves, e) The R-band light curve 
shows hourly average fluxes of the source, the error bars are 
smaller than the symbols in most cases, f) The optical polar- 
ization (filled circles, left axis) and polarization position angle 
(triangles, right axis) are hourly averages. 



tral index as in 2011). The simultaneous Fermi-LAT spectrum 
was calculated for the whole interval from 2010 January to June. 
There was no simultaneous X-ray observation, while for the op- 
tical we use the average (host galaxy subtracted) flux from nights 
when MAGIC was also observing. This "low state SED" is pre- 
sented for illustrative purposes only but was not modeled, since 
both the synchrotron and IC peaks are poorly constrained. 

The SED of 2011 shows two peaks, with the synchrotron 
peak frequency slightly above the optical band, as found for 
many other VHE y-ray emitting BL Lac objects. The X-ray spec- 
tral index is also typical for a BL Lac source. The second peak 
seems to be located between the Fermi-LAT and MAGIC points 
(~1 GeV) as for many of the VHE y-ray emitting BL Lacs. The 
location s of the synchrotron and IC peaks agree with values de- 
rived in lAbdo et alj d2010al) for this source, but the synchrotron 
peak lum inosity was slightly h igher than in the previous obser- 
vation bv lGiommi et alj d2012l) . 

The emission characteristics of BL Lac objects is gen- 
erally well reproduced by the one-zone leptonic model, in 
which a population of relativistic electrons inside a region 
moving down the jet emit throu gh synchrotron and syn- 
chrotron self-Compton mechanisms dBloom & Marschen, 119961 : 
iTavecchio et al.lfl998h . The spectral energy distribution in 201 1 
w as modeled with the on e -zone leptonic model fully described 
in iMaraschi & Tavecchiol d2003l) . The emission region was as- 
sumed to be spherical, with radius R, filled with a tangled mag- 
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Fig. 8. Spectral energy distribution of 2011 January-February 
data (red symbols) modeled with the one-zone SSC model of 
iMaraschi & Tavecchiol d2003l) . From high to low energies: the 
deabsorbed MAGIC spectra (asterisk, see text), the Fermi-LAT 
data (filled squares), Swift-XRT and Swift-UVOT data (trian- 
gles: red for MJD 55569, blue for MJD 55565) and simultane- 
ous KVA data (filled circle, host galaxy subtracted, see text). The 
cyan symbols report the January-June data of Fermi-LAT (data 
points and arrows) and MAGIC (the thick oblique line). Fermi- 
LAT (filled circles) and MAGIC (thick line) data. The green 
open squares are archival data. The dashed line is the model fit 
using the extreme Doppler factor 6 = 60, while the solid line is 
the model fit with high y,„,„ and the long dashed line reports the 
model parameters that produces with smallest x 2 (see text and 
Table 2). 



FiglD As expected from the discussion above, we find a large 
Doppler factor, 6 = 60, well above the typical range of Doppler 
factors obtai ned from the modeling of the emission of similar 
sources (e.g. iTavecchio et all 1201 Oh and disagreeing with the 
lower values required by the FRI-BLLac0 unification scheme 
(lUrrv & Padovanill 19951) . However, there is an alternate, viable 
way to reproduce the observed SED using smaller Doppler fac- 
tors; assume a relatively large minimum Lorentz factor of the 
emitting electrons, ymin = 8 x 10 3 . This, along with a steep 
high energy electron energy distribution («2 — 4.85), allows us 
to properly reproduce the narrow synchrotron bump and to lo- 
cate the SSC peak at high enough energies using a moderately 
large boosting, 6 = 30. This solution resembles one discussed 
for the case of BL Lacs showing hard spectra in the soft X-ray 

120051 ITavecchio etafl 120091: 



and TeV band ([K atarzvriski et al 



Kaufmann et 



1 QKat arzv 



Lefa et al 



2011 



). Interestingly, such pa- 



rameters (large ymi n , steep slope) are consistent with the pre- 
diction of some simulations of part i cle ac c eleration by relativis- 
tic sh ocks (e.g. IVirtanen & Vainiol 120031: ISironi & Spitkovskvl 
1201 ll) . For example, for a proton-electron composition, it is ex- 
pected that the electrons are heated when crossing the shock to a 
typical Lorentz factor of r=r re i m p /m e , where m p /m e = 1836 is 
the proton to electron mass ratio and r re i = 2 - 3 is the relative 
Lorentz factor between the upstream and the downstream flows. 
From this F (that is equivalent to our parameter y m i„), electrons 
are subsequently accelerated, forming a non-thermal tail that is 
well approximated by a steep (« = 3.5) power law. 

The goodness of the fit can be judged by eye or by % 2 - 
minimization procedure. For the fits presented above the "eye 
estimate" was used, as for the latter the systematic errors of the 
data from different instruments are in the key role. However, 
we also tested the auto matic ^-minimization procedure of 
iMankuzhivil et all d201 ll) with estimated systematical errors of 
2%, 10% and 40% for optical-X-ray, GeV y-rays and VHE y- 
rays respectively. The y ;m „ is fixed to same value as in our high 
6 model (10 3 ) to allow easier comparison. The resulting param- 
eters are shown in Table 2 and the fit with long-dashed (dark 
green) line in Fig. 8. The minimal^- 2 fit results in lower 5, but in 
a high y max and rather large emission region radius R compared 
to other fits, but still compatible with the day scale variability 
observed in X rays and optical. 



netic field of intensity B and relativistic electrons, emitting syn- 
chrotron and synchrotron self-Compton radiation. The electrons 
follow a smoothed broken power law energy distribution with 
normalization K between y m j n and y max , with slopes n\ and m 
below and above the break at yt,. The relativistic boosting is fully 
accounted for by the Doppler factor 6. We note, however, that 
one-zone models cannot reproduce the spectrum at the lowest 
frequencies, since the emission is self-absorbed below the mil- 
limeter band. It is generally assumed that this part of the SED 
is due to outer regions of the jet that is not important for the 
modeling of the high-energy emission. 

The optical-UV and X-ray data define a narrow synchrotron 
component peaking around 10 15 Hz. At high energies, the SSC 
bump is well constrained by the Fermi-LAT and MAGIC data 
to peak at about 10 GeV. This particular structure of the SED 
is not easy to reproduce. In particular, the relatively wide sep- 
aration between the two peaks inevitably implies a large value 
of the Doppler factor if standard parameters are used for the 
electron energy distribution ( e.g. |Georganopoulos & Kazanasl 
120031: ITavecchio & Ghiselli"n"iLl2008l) . Our best attempt to repro- 
duce the data in the standard framework provides the param- 
eters given in Table 2 and is displayed as the dashed line of 



5. Summary and Conclusions 

In this paper the first detection of VHE y-rays from 
1ES 1215+303, resulting from MAGIC observations triggered 
by an optical outburst of the source in 2011 January, has been 
reported. In those data, the source is clearly detected at a 9.4 cr 
significance level. Also simultaneous multi-wavelength data are 
presented from radio to HE y-rays and compared to results from 
earlier MAGIC observations in 2010, when the source was in a 
lower optical state. The VHE y-ray flux in 201 1 was higher com- 
pared to 2010, suggesting that the activity in these two bands 
is connected. This conclusion is further supported by the fact 
that 1ES 1215+303 is already the fifth discovery at VHE y-rays 
achieved after the MAGIC observations were triggered by an op- 
tical outburst. 

Our collected multi-wavelength data set is the most extensive 
energy coverage for 1ES 1215+303 to date. The optical- VHE 
y-ray outburst seems to have been accompanied by an X-ray 
outburst, while in the Fermi-LAT band the flux increased only 
marginally. The optical photo-polarimetric data suggests that the 

6 Fanaroff-Riley I radio galaxies (FR I) 
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Table 2. Input model parameters for the three models shown in Fig|8] The following quantities are reported: the minimum, break, 
and maximum Lorentz factors and the low and high energy slope of the electron energy distribution, the magnetic field intensity, 
the electron density, the radius of the emitting region and its Doppler factor. In addition in the last column we report the^ 2 /d.o.f 
assuming 2%, 10% and 40% systematical errors for optical-X-ray, GeV y-rays and VHE y-rays respectively. 



high state could be caused by a shock traveling down the jet that 
collides with a standing shock with a differently oriented mag- 
netic field. The X-ray and VHE y-ray high states could then also 
originate from this collision. 

The SED of 1ES 1215+303 in 2011 January was modeled 
using a one-zone SSC model since it provides a good descrip- 
tion of the SED of many VHE y-ray emitting BL Lac objects. 
However, for 1ES 1215+303 the synchrotron and IC peaks are 
narrow, the separation between the two peaks is wide, and a sim- 
ple one-zone SSC model with typical parameters failed to repro- 
duce the observed SED. To fit the SED, a high Doppler factor 
or a narrow electron energy distribution is required. While high 
Doppler factors are disfavored by the unified models, the high 
y„„„ value could be a viable solution in the light of simulations 
modeling the acceleration of electrons in a relativistic shock in 
a proton-electron jet. This should be further investigated, e.g. 
using the fully self-consistent SSC model with particle accelera- 
tion due to shock and stochastic acce leration bv lWeidinger et al] 
(l2010h : IWeidinger & Spaniel (1201 Oh . 

Given the rather extreme conditions needed for the one- 
zone model, the presence of a ve l ocity structure i n the jet 
dGeorganopoulos & KazanasL l2003t iGhisellini etail l2005h is 
also possible and should be tested for modeling the SED. The 
narrow synchrotron and IC peaks should well constrain the 
model. Additionally, the more complex emission scenario sug- 
gested by the photo-polarimetric behavior o f 1ES 12 1 5+30 3 
should be tested using, for instance, a model o flMarsched ( l20Tlh . 
who investigated the emission from a turbulent ambient jet 
plasma that passes through either the standing or moving shocks 
in the jet. Further observations of this exceptional VHE y-ray 
emitting BL Lac are strongly encouraged. 
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